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ABSTRACT 

We investigate the limitations of thermonuclear X-ray bursts as a distance indicator for the weakly- 
magnetized accreting neutron star 4U 1728—34. We measured the unabsorbed peak flux of 81 bursts in 
public data from the Rossi X-Ray Timing Explorer. The distribution of peak fluxes was bimodal: 66 
bursts exhibited photospheric radius expansion (presumably reaching the local Eddington limit) and were 
distributed about a mean bolometric flux of 9.2 x 10~ 8 ergcm~ 2 s _1 , while the remaining (non-radius 
expansion) bursts reached 4.5 x 1CU 8 ergcm _2 s _1 , on average. The peak fluxes of the radius-expansion 
bursts were not constant, exhibiting a standard deviation of 9.4% and a total variation of 46%. These 
bursts showed significant correlations between their peak flux and the X-ray colors of the persistent 
emission immediately prior to the burst. We also found evidence for quasi-periodic variation of the 
peak fluxes of radius-expansion bursts, with a time scale of ~ 40 d. The persistent flux observed with 
RXTE/ ASM over 5.8 yr exhibited quasi-periodic variability on a similar time scale. We suggest that these 
variations may have a common origin in reflection from a warped accretion disk. Once the systematic 
variation of the peak burst fluxes is subtracted, the residual scatter is only ~ 3%, roughly consistent with 
■^j- ' the measurement uncertainties. The narrowness of this distribution strongly suggests that i) the radiation 

\Q , from the neutron star atmosphere during radius-expansion episodes is nearly spherically symmetric, and 

ii) the radius-expansion bursts reach a common peak flux which may be interpreted as a standard candle 
intensity. Adopting the minimum peak flux for the radius-expansion bursts as the Eddington flux limit, 
we derive a distance for the source of 4.4-4.8 kpc (assuming i?NS = 10 km), with the uncertainty arising 
from the probable range of the neutron star mass Mns = 1.4-2 M . 

Subject headings: stars: neutron — X-rays: bursts — nuclear reactions — equation of state — stars: 
individual (4U 1728-34) — X-rays: individual (4U 1728-34) 
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C$ • changes in the X-ray intensity of accreting neutron stars in 

low-mass X-ray binary (LMXB) systems, with rise times where Mm ig thfi m&gg of thg neutron T ig the rf _ 

between $1-10 s and decay times between ~ 10 - 100 s fective temperature of the atm0 sphere, a T is a coefficient 
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et al. 2001) and are caused by unstable nuclear burning of . t ,• mi • ,i , , I'nn 
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accreted matter on the neutron-star surface. r • u j ^ j j-i, i. t ± c 

Tr . . . for cosmic abundances), and the parameter f accounts for 

If the thermonuclear energy during a burst is released pogsible anisotropy of the burst emission. The final fac- 

sufficiently rapidly, the flux through the neutron star at- tQr ^ parcntheses repreS ents the gravitational redshift due 

mosphere may reach the local Eddington limit, at which tQ thfi com t nature of the neutron gta and also de _ 
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ergy of the X-ray photosphere, which is lifted above the ity dependg Qn the ratio JVf NS J i? NS , measurements of the 

neutron star surface while the^ emerging luminosity (mea- peak flux of radius _ expans ion bursts allows in principle the 

sured locally) remains approximately constant and equal me asurement of these fundamental properties (e.g. Damen 

to the Eddington limit. These are the so-called radius- gt &1 lggo Smale 20Q1 Kuulkerg et aL 2002b) , Qn the 

expansion or Eddington- limited bursts. , , u j c j. ui i 

• ■ • , n i i- ^ other hand, because the masses and radii of stable neutron 

For spherically symmetric emission, _ the Eddington hi- gtarg predicted b iven equation of state span a nar _ 

mmosity measured by an observer at infinity is given by row r£mge Qf yalueg ^ e g > Lattimer & Prakash 20 01), 

(Lewm et al. 1993) Eddington-limited X-ray bursts can be used as distance 
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^Edd,oo — Tj- — ^ I 1 I The validity of the physical picture of Eddington-limited 
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bursts discussed above can be verified observationally in 
two ways. First, the peak fluxes of such bursts from 
each individual source should be the same. Second, the 
peak luminosities inferred for sources with independent 
distance measurements, such as those in globular clusters, 
should correspond to the Eddington limit for a neutron 
star. Since the original discovery of Eddington-limited 
bursts, a number of authors have addressed these ques- 
tions. In early observations, a significant number of radius 
expansion bursts had been detected from three sources 
and their peak fluxes were found to be similar to within 
~ 20% (4U 1820-30: Vaccaet al. 1986; Damenet al. 1990; 
4U 1636-536: Damcn ct al. 1990; 4U 1728-34: Basinska 
et al. 1984). In particular, the peak luminosities of radius- 
expansion bursts observed from 4U 1820—30, which resides 
in the gobular cluster NGC 6624 were found to be compa- 
rable to the Eddington limit for a neutron star. Although 
successful in providing support to the model of Eddington- 
limited bursts, these early studies were limited to a small 
number of sources and suffered from the statistical uncer- 
tainties inherent to fitting spectral models to low signal-to- 
noise data. In recent years, observations with BeppoSAX 
and the Rossi X-ray Timing Explorer (RXTE) have re- 
vealed a large number of Eddington-limited bursts from 
several sources, which have been studied in great detail. 
Using BeppoSAX and RXTE data, Kuulkers et al. (2002a) 
recently studied the Eddington-limited bursts of 12 globu- 
lar cluster sources with well known distances and showed 
that, with one exception, their peak fluxes were constant 
to within ~ 15% and were comparable to the Eddington 
limit for neutron stars with H-poor atmospheres. 

In this series of articles, we use all the publicly avail- 
able data obtained with RXTE to date in order to ob- 
servationally test the hypothesis that Eddington-limited 
bursts can be used as distance indicators for neutron-star 
LMXBs. In the present study we quantify, and examine 
the causes of, systematic variations in the peak burst flux 
from the source with the greatest number of bursts de- 
tected by RXTE, 4U 1728-34. 



2. RXTE OBSERVATIONS OF 4U 1728-34 

The X-ray source 4U 1728-34 (GX 354+0; I = 354?3, 
b = — 0?15) was first resolved by Uhuru scans of the Galac- 
tic center region (Forman et al. 1976). Thermonuclear X- 
ray bursts from 4U 1728—34 were discovered during SAS-3 
observations (Lewin et al. 1976; Hoffman et al. 1976). The 
bursting behavior was subsequently studied in detail us- 
ing extensive observations by SAS-3, which accumulated 
96 bursts in total. From these data Basinska et al. (1984) 
showed evidence for a narrow distribution of peak burst 
fluxes, as well as a correlation between the peak flux and 
the burst fluence. The distance to the source has pre- 
viously been estimated from measurements of the peak 
burst fluxes as between 4.2-6.4 kpc (van Paradijs 1978; 
Basinska ct al. 1984; Kaminker et al. 1989). The esti- 
mated extinction to the source is Ay « 14; only a precise 
position following from the detection of a radio counter- 
part allowed identification with a K = 15 infrared source 
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Fig. 1.— RXTE/ ASM measurements of 4U 1728-34. The top 
right panel shows the Lomb-normalized pcriodogram computed over 
the entire ASM history. The top left panel shows the corresponding 
dynamical Lomb pcriodogram, calculated from 730-d subsets of the 
ASM data at 9-d intervals. The lower panel shows the 5-d averaged 
count rates, computed from the initial 1-d averages excluding those 
points where the measured error > 0.5 counts -1 . Error bars, where 
shown, indicate the 1<t uncertainties. 



(Marti et al. 1998). No independent distance measurement 
is available. Long-term Ariel-5 measurements, as well as 
the first two years of monitoring by the All-Sky Monitor 
(ASM) on board RXTE, suggest the presence of a long- 
term quasi-periodicity of 63 or 72 d, respectively (Kong 
et al. 1998). 

RXTE/PCA observations of the source in 1996 led to 
the discovery of nearly coherent millisecond oscillations 
during the X-ray bursts (Strohmayer et al. 1996). Similar 
oscillations were subsequently observed in 9 other sources 
(van der Klis 2000; Wijnands et al. 2001; Galloway et al. 
2001; Kaaret et al. 2002). A substantial archive (1140 ks) 
of public RXTE/PCA data from 4U 1728-34 has accumu- 
lated throughout these observations, dating from shortly 
after the launch of the satellite on 1995 December 30. Sub- 
sets of the bursts observed during the PCA observations 
have been studied by van Straaten et al. (2001) and Franco 
(2001), with particular attention to the relationship be- 
tween the appearance of burst oscillations and the mass 
accretion rate. A significant fraction of the bursts observed 
by RXTE show evidence for photospheric radius expan- 
sion, exhibiting the characteristic temporary increase in 
the apparent blackbody radius, coincident with a decrease 
in the color temperature in the early stages of the burst. 
Muno ct al. (2001) used this dataset to find a correlation 
between the frequency of the burst oscillations and the 
preferential appearance in radius expansion or non-radius 
expansion bursts in nine sources. 

We have obtained all the available public RXTE 
data to date from the High-Energy Astrophysics 
Science Archive Research Center (HE AS ARC; 
http://heasarc.gsfc.nasa.gov). This study is part of a 
larger effort involving more than 70 known bursters; the 
burst detection and analysis procedures, as well as the 
resulting burst database are described in more detail in 
Galloway et al. (2003b, in preparation). Several of the 
bursts from 4U 1728—34 we analyzed have also been stud- 
ied in detail by Strohmayer et al. (1996, 1997, 1998). 
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2.1. ASM 

The All-Sky Monitor (ASM) onboard RXTE consists of 
three Scanning Shadow Cameras (SSCs) sensitive to pho- 
tons in the energy range 2-10 keV, mounted on a rotating 
platform (Levine et al. 1996). ASM observations are per- 
formed as sequences of dwells lasting 90 s, after which the 
platform holding the SSCs is rotated by 6 degrees. Most of 
the sky is viewed once every few hours. The data from each 
SSC from each dwell are averaged to obtain the daily inten- 
sities of known sources in the field of view. The long-term 
2-10 keV RXTE/ ASM flux history of the source is shown 
in Fig. 1. The observed photon flux exhibited variations 
of <~ 5 count s -1 on time scales of ~ 10 d, superimposed 
on a long-term trend of decreasing mean intensity. 

We computed the Lomb-normalised periodogram (Press 
& Rybicki 1989) of the full dataset, shown also in Fig. 1. 
We found significant evidence for a periodicity with 
Pasm = 38.6 d, with a Lomb-normalized power of 31.9 
as well as secondary peaks at 63.7 and 67.5 d. The latter 
values are close to those found from earlier Ariel-5 mea- 
surements, as well as a 2-yr subset of the RXTE/ ASM data 
by Kong et al. (1998). The rms amplitude of the signal 
at Pasm was ~ 8 — 9%. The dynamical Lomb-normalized 
periodogram shows that this periodicity varied in strength 
over the entire ASM history (Fig. 1). The « 39 d periodic- 
ity appeared strongly only in the first half of the measure- 
ment interval, when the source reached its peak long-term 
intensity. After MJD 51400, the dominant periodicity ap- 
peared to shift to around 30 d. A broad peak around 65 d 
was also present at most times. 

2.2. PCA 

The Proportional Counter Array (PCA; Jahoda et al. 
1996) consists of five identical gas-filled proportional 
counter units (PCUs) with a total effective area of « 
6000 cm 2 , and sensitivity to X-ray photons in the 2-60 keV 
range. We initially scanned 1 s binned lightcurves over the 
full PCA energy range (created from "Standard-1" mode 
data) in order to locate X-ray bursts from 4U 1728—34. 
In all the PCA observations, data were also collected in 
user-defined modes offering higher time resolution, com- 
pared to the standard data modes. The PCA records the 
arrival time (1 fis resolution) and energy (256 channel res- 
olution) of every unrejected photon; data were generally 
binned on somewhat lower time and spectral resolution in 
order to meet telemetry limits. Once an X-ray burst was 
detected, we extracted 2-60 keV PCA spectra within in- 
tervals of 0.25-1 s covering the burst. A response matrix 
was generated separately for each burst using pcarsp ver- 
sion 8.0 (part of lheasoft release 5.2, 2002 June 25) in 
order to take into account the known gain variations over 
the life of the instrument. The gain was manually re-set by 
the instrument team on 3 occasions (1996 March 21, 1996 
April 15 and 1999 March 22); an additional gain epoch (5) 
began on 2000 May 13 with the loss of the propane layer 
in PCU #0. In addition to these abrupt changes, a more 
gradual variation in the instrumental response is known to 
occur, due to a number of factors. 

We then analyzed these data using an approach that is 
often used in X-ray burst analysis (e.g. Kuulkcrs et al. 
2002b, although see also van Paradijs & Lewin 1986). For 
the background to the burst spectra we used a spectrum 
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Fig. 2. — The peak flux -F pe ak (excluding the pre-burst persistent 
emission) plotted as a function of the integrated burst fluence E\, 
of 15 non-radius expansion (cross) and 66 radius expansion (dia- 
mond) bursts from 4U 1728—34. The result of a linear fit to F pea k 
as a function of for the radius expansion bursts is plotted as 
a dashed line. Error bars indicate the la uncertiainties. Note the 
broken x-axis. 

extracted from a (typically) 16 s interval prior to the burst. 
Each time-resolved background-subtracted spectrum dur- 
ing the bursts was fit with a blackbody model multiplied 
by a low-energy cutoff representing interstellar absorption 
with fixed abundances. The initial spectral fitting was per- 
formed with the absorption column density nn free to vary; 
the resulting fit values typically exhibited very large scat- 
ter, particularly towards the end of the burst when the flux 
was low. Thus, for the final analysis we re-fit each spectra 
with nn fixed at the weighted mean value measured over 
the entire burst. The unabsorbed bolometric flux .Fboi,i a t 
each timestep was then estimated according to 

= 1.0763 x 10- 11 T^K^ erg cm" 2 s" 1 (2) 

where i?NS is the neutron star radius, d is the distance 
to the source, Tbb is the blackbody (color) temperature 
in units of keV, and A'bb is the normalisation of the 
blackbody component as returned by the fitting program 
(xspec version 11). We also estimated the burst fluence 
Eb by integrating the measured Fb i,» over the burst dura- 
tion. We discuss the possible consequences of the bolomet- 
ric correction implicit in equation (2) in section §3.1. Ku- 
ulkcrs et al. (2002b) and other authors have noted a « 20% 
systematic flux offset in RXTE measurements compared 
to other instruments. It may be argued that the absolute 
flux calibration of these instruments is no better than that 
of RXTE: in any case, the RXTE data still offer substan- 
tially better signal-to-noise and hence greater precision of 
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Fig. 3. — The peak fluxes -F pca k (excluding the pre-burst persis- 
tent emission) of 81 X-ray bursts from 4U 1728—34 as a function of 
burst number, which increases monotonically with time. The hori- 
zontal solid line shows the mean peak flux of the radius expansion 
bursts (diamonds), while the dashed lines show the lcr limits. Error 
bars indicate the lcr uncertainties on each measurement. 
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Fig. 4. — The Lomb- normalized periodogram of the time vari- 
ation of the peak fluxes Fp ea k re from 50 radius-expansion bursts 
occurring before MJD 51500 ("of 66 radius expansion bursts in to- 
tal). Note the indications for excess power between 30-60 d; the 
most significant peak is at 38 d, with a Lomb power of 16.3. 

flux measurements. Throughout this paper we quote the 
unadjusted RXTE peak fluxes, but for distance estimates 
(see §3.4) we consider the possible effects of this systematic 
offset. 



To measure the (pre-burst) persistent flux we estimated 
the instrumental background using pcabackest version 
3.0, and fit the background-subtracted pre-burst spectrum 
with an absorbed blackbody and power law model. The 
typical reduced-% 2 (xV) was 1-2- Gain-corrected X-ray col- 
ors were also calculated from the pre-burst spectrum. We 
determined correction factors for the count rates in the 
various energy bands by comparison of Crab spectra over 
the lifetime of the instrument (see Muno et al. 2002b, for 
more details of the correction method) . The soft color was 
calculated as the ratio between the background-subtracted 
gain-corrected count rates in the 3.6-5.0 keV and 2.2- 
3.6 keV bands, and hard color as the ratio between counts 
in 8.6-18 keV and 5.0-8.6 keV bands. 

2.3. Peak burst fluxes 

We detected 81 X-ray bursts observed by RXTE from 
4U 1728-34 between 1996 February 15 and 2001 Novem- 
ber 15, with peak burst fluxes in the range (0.3-1.2) x 
10~ 7 ergcm _2 s _1 . We also detected four much fainter 
bursts, which peaked at w 5 x 10~ 9 ergcm _2 s _1 ; two on 
1996 May 3, and one each on 2001 May 27 and May 29. 
These burst times closely followed active periods of the 
nearby source MXB 1730—33 (the "Rapid Burster"), just 
0.53° away, compared to field of view of RXTE of 1°. Since 
the next brightest burst from 4U 1728—34 was a factor of 6 
brighter, we conclude that the four very faint bursts actu- 
ally originated with the Rapid Burster, and exclude them 
from our analysis. 

Most of the X-ray bursts from 4U 1728-34 exhibited 
some degree of radius expansion. As a working definition, 
we considered that the local Eddington limit had been 
reached when (i) the blackbody normalization K^, which 
is directly related to the surface area of the emitting region, 
reached a local maximum close to the time of peak flux; 
(ii) lower values of the normalization were measured 
following the maximum, with the decrease significant to 
4cr or more; and (iii) there was evidence of a significant 
(local) decrease in the fitted temperature Tbb at the same 
time as the increase in the normalization ifbb- For 8 of 
the bursts there was weak evidence of radius expansion, 
at only the 2-3cr level. However, the distribution of peak 
fluxes in this class of bursts was consistent with that of the 
bursts exhibiting more significant radius expansion, and so 
we treated them as one population. 

The peak bolometric fluxes -Fpeak^ rnax(i 7 boi.i) were 
bimodally distributed, with the radius-expansion bursts 
around a mean of 9.2 x 10~ 8 ergcm _2 s _1 , and the non- 
radius expansion bursts around 4.5 x 10~ 8 crgcm _2 s _1 
(Fig. 2). The peak fluxes of the radius expansion bursts 
were roughly proportional to the fluence, up to an asymp- 
totic level of w 1.2 x 10~ 7 crgcm~ 2 s _1 (Spearman's rank 
correlation p = 0.78, significant at approximately the 
6cr confidence level). The burst fluence generally ranged 
from (0.04-0.9) x 10~ 6 ergcm~ 2 . The non-radius ex- 
pansion bursts had fluences below 0.5 x 10~ 6 crgcm~ 2 , 
while the majority of the radius expansion bursts had flu- 
ences between (0.4-0.8) x 10~ 6 crgcm~ 2 . The brightest 
radius-expansion burst (#72, on 2001 February 9 03:01.40 
UT) also had the largest fluence. While its peak flux at 
(1.21 ± 0.05) x 10~ 7 ergcm~ 2 s" 1 was around 30% larger 
than the mean for the remaining radius expansion bursts, 
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the fluence was almost 150% greater than the mean at 
(1.409 ±0.016) x 10~ 6 erg cm" 2 . 

In the radius expansion bursts, the peak flux was gener- 
ally reached during the contraction stage following the ra- 
dius maximum. This is contrary to basic theory of Edding- 
ton limited bursts, which predicts that the observed bolo- 
metric flux should be constant throughout the radius ex- 
pansion and contraction (once the flux is corrected for the 
effect of gravitational redshift, which is naturally variable 
due to the changing elevation of the photosphere above 
the neutron star surface). Furthermore, as noted by van 
Straaten et al. (2001), some bursts exhibited unusual vari- 
ation of the fitted radius as the burst evolved. In about 
half the bursts the radius expansion and subsequent con- 
traction were highly significant, but following the radius 
minimum (which is usually assumed to be the time when 
the expanding material "touches down" on the neutron 
star surface) the apparent radius increased again, up to a 
level which in some cases was as high as the initial radius 
maximum. Since the initial rise and fall in the fitted ra- 
dius generally exceeded our criteria for classification as a 
radius expansion burst, we included these unusual cases 
with the other "normal" radius expansion bursts. How- 
ever, there are clearly additional factors influencing the 
observed spectra throughout the burst (possibly involving 
changes in the photospheric composition, perhaps related 
to ejection of an hydrogen-rich envelope; Sugimoto et al. 
1984). While these additional effects may also help to de- 
termine the peak flux of radius expansion bursts, in gen- 
eral their greatest influence appears to be on the evolution 
later in the burst, after the initial 5 s or so when the radius 
expansion takes place. 

2.4. Peak flux variation in radius- expansion X-ray bursts 

The peak fluxes of the radius-expansion bursts from 
4U 1728-34 show a significant ( X 2 = 1760 for 60 de- 
grees of freedom) deviation from a constant value (Fig. 
3). The standard deviation of the peak fluxes was 8.7 x 
10~ 9 ergcm _2 s _1 , corresponding to a fractional rms of 
9.4%; the net variation was 46%. 

While the bursts were clustered extremely irregularly in 
time (primarily due to the irregular scheduling of pointed 
observations) the peak fluxes of radius expansion bursts 
that were close in time gave a strong suggestion of system- 
atic evolution, on time scales of a few tens of days. To illus- 
trate this evolution, we calculated a Lomb-normaliscd pe- 
riodogram of the peak flux of the radius-expansion bursts. 
Like the ASM periodicity, the time scale of the variation 
in the peak fluxes did not appear to be consistent over 
the full set of bursts we measured from RXTE observa- 
tions. Thus, we selected only those bursts before MJD 
51500 to calculate the periodogram (note that this is also 
when the ~ 40-d periodicity in the ASM lightcurve be- 
came much weaker; Fig. 1). The periodogram shows ev- 
idence for excess power between 30-60 d; the most sig- 
nificant peak was at 38 d, with a Lomb power of 16.3 
(estimated significance > 5er from Monte-Carlo simula- 
tions; Fig. 4). We fit the peak fluxes folded on the 
38 d period with a sine curve, with a resulting fractional 
rms amplitude of 7.7%. The next most significant peak 
in the periodogram was at 26 d, with power 9.0 (3.1er). 
Additional evidence that the changes in the peak burst 
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Fig. 5. — Peak flux -F pe ak,RE °f radius-expansion thermonuclear 
X-ray bursts from 4U 1728—34, as a function of the gain-corrected 
soft (3.5-5.0 keV/2.2-3.6 kcV) and hard (8.6-18 keV/5.0-8.6 keV) 
X-ray colors of the persistent flux, prior to each burst. Error bars 
represent the estimated la uncertainties. 



flux was linked to the long-term source evolution comes 
from significant correlations measured between -F pea k,RE 
and both the hard and soft gain-corrected colors (Fig. 5). 
The (Spearman's) rank correlation between the soft color 
and -Fpeak.RE was p = 0.62, with a significance of 6.8 x 10~ 6 
(equivalent to 4.5cr); for the hard color, p = 0.84, with a 
significance of 3.4 x 10~ 15 (equivalent to 6.1<r). 

While the 38 d periodicity found in the radius-expansion 
burst peak fluxes was corroborated by a significant de- 
tection at a similar period in the ASM flux history, the 
26 d period was not. Thus, we consider the evidence 
for the presence of the latter signal to be weak. The ap- 
parent quasi-periodic behavior of the ASM flux modula- 
tion, coupled with the similarity in the time scales and 



6 



Galloway et al. 




S3 50120 50140 51110 51130 

OS 

I 




-10 



51190 51210 51230 51250 
Time (MJD) 

Fig. 6. — The peak burst fluxes -F pea k,R,E of selected subsets of 
radius-expansion bursts from 4U 1728—34, with the best-fit second- 
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ror bars show the la uncertainties. 
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Fig. 7. — The combined distribution of residuals from the fits 
shown in Fig. 6, with the expected Gaussian distribution for the 
measured a = 2.8% overplotted. 

modulation amplitudes of the ASM flux and the radius- 
expansion burst peak fluxes, suggest that these two phe- 
nomena share a common cause. The uneven time sampling 
of -Fpeai^RE and the time-dependence of the ASM period- 
icity (Fig. 1) did not allow us to test this hypothesis for 
the entire set of radius-expansion bursts. However, we 



have identified four intervals, in which multiple radius ex- 
pansion bursts were detected over time spans of no more 
than « 30 d. The resulting 4 subsamples contain a to- 
tal of 41 of the 66 radius-expansion bursts measured from 
4U 1728—34. For each of these intervals we attempted 
to remove the long-term variation by fitting the overall 
-Fpcak.RE time-dependence using a second order polynomial 
(Fig. 6). The combined distribution of the residuals from 
each of the polynomial fits is shown in Fig. 7. The dis- 
tribution is consistent with Gaussian noise with a = 3.2% 
of the mean peak flux, compared to the averaged la flux 
measurement errors of 2.1%. 

3. DISCUSSION 

We have studied 81 thermonuclear X-ray bursts from 
4U 1728—34, which includes 66 exhibiting evidence 
for photosphcric radius expansion, as observed by the 
RXTE/PCA. We have shown that the radius expansion 
bursts exhibit a significant variation in their peak fluxes 
-Fpcak.RE- The fractional standard deviation was ~ 9.4%, 
while the total variation was ~46%. This is significantly 
larger than the formal measurement errors, which were 
typically ~ 2%. This result appears inconsistent with the 
simple picture of radius-expansion bursts (e.g. Lewin et al. 
1993), according to which the peak flux should be nearly 
constant and equal to the Eddington critical flux. How- 
ever, we must also consider the possibility of a variety of 
systematic effects arising as a consequence of our analysis 
method, which may contribute (or give rise) to the mea- 
sured variation. 

3.1. Possible systematic effects 

While subtraction of the pre-burst emission as back- 
ground has dubious theoretical justification, and may have 
originally been adopted for convenience, the method has 
been shown in several analyses to be relatively robust 
(e.g. van Paradijs & Lewin 1986; Kuulkers et al. 2002b). 
An implicit assumption is that the persistent emission re- 
mains unchanged throught the burst. This may not be 
true, especially for the very energetic, radius-expansion 
bursts that may disrupt the inner accretion flow. As a 
result, variations in the true persistent flux F pcr (t) dur- 
ing the burst may contribute to variation in the measured 
bolometric burst flux. The persistent flux prior to the ra- 
dius expansion bursts in 4U 1728—34 was typically around 
3.7 x 10~ 9 ergcm~ 2 s -1 , but for some bursts was as high 
as 5 x 10 -9 ergcm~ 2 s _1 . The net variation over all the 
bursts was only 3.2 x 10~ 9 ergcm _2 s _1 , which is insuf- 
ficient to account for the (4.2) x 10~ 8 ergcm _2 s _1 net 
variation in the peak burst fluxes. Even if the persistent 
flux was quenched completely during the burst, this would 
still only give a variation of at most 5 x 10~ 9 erg cm~ 2 s _1 . 
Thus, variations in the peristcnt flux cannot completely ac- 
count for the observed variation in .Fpeak.RE, but may con- 
tribute to some extent. If this is the case, we may expect 
a correlation between the F pea k,RE and the pre-burst per- 
sistent flux. In fact, -Fpoak.RE and the persistent flux were 
weakly anticorrelated, with Spearman's rank correlation 
p = —0.32 (estimated signifcance 1.5 x 10~ 2 , equivalent to 
2.4<j). If anything, the persistent flux variation serves to 
suppress slightly the true variation in -F p0 ak,RE- 
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The bolometric correction (equation 2) typically adds 
~ 7% to the peak 2.5-20 keV flux of the radius expansion 
bursts as measured by the PCA. This correction may in- 
fluence our result in one of two ways. On the one hand, 
systematically biased bolometric corrections may give rise 
to a variation in flux which is not present in the observed 
source fluxes (restricted to the PCA passband). We can 
easily rule out this possibility, since the flux integrated 
just over the PCA passband exhibits identical (fractional) 
variation as the inferred bolometric flux. Alternatively, 
deviations in the spectrum from a perfect blackbody may 
give rise to (unobserved) variations in the flux outside the 
PCA passband that contribute to the true bolometric flux 
variation. In the latter case these unobserved bolometric 
flux variations may (for example) partially or completely 
compensate for the variation we see in the PCA passband, 
so that the bolometric flux variation we infer is exagger- 
ated, or erroneous. The second effect is more difficult to 
discount, since (obviously) the spectral variations outside 
the PCA passband are not measurable with the present 
data. However, the ~ 7% typically added to the flux in 
the PCA passband by the bolometric correction is signif- 
icantly smaller than the ~46% observed variation in the 
peak fluxes. In order to compensate for the flux variation 
we observe, this contribution would then need to fluctuate 
by a factor of ~ 6, which seems unlikely. 

Thus, we conclude that the observed variation was not 
an artifact of any aspect of our analysis method and, 
hence, must be genuine. 

3.2. Origin of the peak flux variation 

The large number of radius-expansion bursts observed 
from 4U 1728-34 with RXTE/PCA along with the long- 
term flux history accumulated by the ASM provides for 
the first time a plausible cause for the observed broad 
distribution of peak burst fluxes. The fact that (i) the 
peak burst fluxes were correlated with the X-ray colors of 
the persistent emission, (ii) they varied in a quasi-periodic 
manner, and (iii) the timescale and fractional amplitude of 
the variability were similar to those of the persistent emis- 
sion, strongly suggest that the same phenomenon causes 
the variability in both the persistent and burst fluxes. 

Since the variability of the persistent emission is not 
coherent, it is unlikely to be due to orbital modulation. 
As noted by Kong et al. (1998), its timescale (whether 
~ 30 or 60-70 d) is much longer than would normally 
be expected for the orbital period of a Roche lobe-filling 
low-mass X-ray binary. Such "super-orbital" periodicities 
are observed in several similar sources (e.g. White et al. 
1995), and are generally attributed to variations in the 
accretion geometry, possibly caused by the precession of 
a warped accretion disk about the neutron star. If the 
persistent emission is modulated at this timescale because 
of a slowly evolving warp in the accretion disk that is re- 
flecting a small fraction of the X-ray luminosity of the 
central object to the observer, then the peak burst fluxes 
would also be modulated at the same timescale and with 
a similar amplitude. Our analysis of the radius-expansion 
bursts of 4U 1728—34 strongly suggest that this is the case. 
Further support is provided by the highly significant cor- 
relation observed between the peak flux and the fluence 
for the radius-expansion bursts. If there was no additional 



contribution to the burst flux from disk reflection, theory 
predicts that the peak flux would be independent of the 
fluence. 

It remains to be established whether, given the condi- 
tions in the 4U 1728—34 system, the accretion disk can 
become warped, undergo precession at approximately the 
measured (quasi-) periodicities, and give rise to the ob- 
served degree of modulation of the persistent and peak 
radius-expansion burst fluxes. The disk warp may arise 
(for example) from non-axisymmetric radiation pressure 
forces (e.g. Pringle 1996; Maloney et al. 1996; Maloney 
& Begelman 1997). The conditions required for the initial 
warping, and steady precession thereafter, depend primar- 
ily upon the orbital separation and the efficiency of accre- 
tion (Ogilvie & Dubus 2001). The orbital parameters are 
presently unknown for 4U 1728—34, but are in principle 
measurable; the accretion efficiency is much more diffi- 
cult to measure for this, or any other, LMXB. Thus, our 
present level of knowledge is not sufficient to reject the 
hypothesis that a precessing, warped disk (whether aris- 
ing by radiation instabilities, or some other mechanism) is 
present in 4U 1728—34. While the question of whether a 
disk warp can give rise to the observed modulation is com- 
paratively more straightforward, we are still limited by the 
lack of measured system parameters, in particular the in- 
clination. For a flat accretion disk Lapidus & Sunyaev 
(1985) calculated an anisotropy factor of 2.8, depending 
upon the inclination angle. The precessing of a warped 
disk is likely to affect the proportion of reprocessed radia- 
tion observed during the burst in the same way that vary- 
ing the inclination would. The derived anisotropy factor 
is more than sufficient to explain the observed modulation 
in the peak flux of radius-expansion bursts. 

Because of these uncertainties, we can most likely adopt 
a relatively wide range of parameters (disk warping angle, 
disk albedo) which will give rise to a modulation of at 
least the level measured in 4U 1728—34; thus, such an ap- 
proach would also have no ability to rule out warped disk 
precession as a mechanism for the X-ray flux modulation. 
We note that in the archetypical precessing warped disk 
system Her X-l, periodic obscuration of the neutron star 
by the disk gives rise to a modulation of the persistent X- 
ray flux of essentially 100% (e.g Scott et al. 2000), much 
larger than the ~ 10% measured for 4U 1728—34. Since it 
exhibits neither X-ray eclipses or dips, 4U 1728—34 must 
have a lower inclination (i < 85°) than Her X-l, making 
obscuration by the disk less likely. For a disk warped to 
the degree inferred for Her X-l (20° at the outer edge), 
the a priori probability for obscuration in 4U 1728—34 is 
~ 25%. Even if the inclination is not sufficiently high to 
permit obscuration, X-ray reflection from the disk, cou- 
pled with the variations in the projected disk area due to 
the precessing warp, may yet be sufficient to give rise to 
the observed modulation. 

3.3. Possible anisotropy of the burst emission 

When the systematic trends in the variation of the peak 
burst fluxes are removed, the residual variation is only 2.8- 
3.2%, which is comparable to the typical measurement un- 
certainty of 2%. This has a very important implication for 
the anisotropy of radius-expansion bursts in 4U 1728—34, 
as described by the parameter £ in equation (1). The small 
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residual scatter of the peak fluxes strongly suggests that 
the intrinsic variation of the peak burst flux is also small, 
~ 1-2%. It seems unlikely that we observe the same face 
of the neutron star at the same orientation during every 
one of these bursts, particularly given the rapid rotation 
inferred from the burst oscillations (364 Hz; Strohmayer 
et al. 1996). Additionally, these same oscillations are al- 
most never observed during the radius expansion episode 
itself, even if they are present earlier or later in the burst 
(Muno et al. 2002a). We conclude that the longitudinal 
dependence of the burst flux during the radius expansion 
episodes is negligible. A latitudinal variation in flux re- 
mains plausible, particularly since the effective gravity is 
smaller at the neutron star equator than at the poles, and 
so we might expect a greater degree of expansion of the 
atmosphere there. However, we observe significant varia- 
tion in the blackbody normalization when the peak burst 
flux is achieved, which suggests that the radius expan- 
sion episodes reach different peak radii. Since we might 
expect the degree of latitudinal anisotropy to vary with 
increasing radius, the effect of such a latitudinal varia- 
tion of flux would be a dependence of the peak flux on 
the blackbody normalisation at the peak, which is not 
observed. Thus, we conclude that the degree of latitu- 
dinal flux anisotropy is most likely also limited by the (in- 
ferred) intrinsic variation of the peak burst fluxes. We 
conclude that the burst emission during the radius expan- 
sion episode is isotropic to within ~ 1-2%. Note that 
it is still possible for the burst emission at the neutron 
star surface to be significantly anisotropic, but that this 
anisotropy is smoothed out through reprocessing in the 
extended atmosphere present during the radius expansion 
episodes. 

3.4. Consequences for distance estimates 

Studies such as this provide a measure of the systematic 
uncertainties of the distance estimates of X-ray sources 
that are based solely on Eddington-limited bursts (see, 
e.g., van Paradijs & White 1995). We note that the stan- 
dard deviation we measure is within the typical peak burst 
flux uncertainty (~ 15%) measured by Kuulkers et al. 
(2002a) for the globular cluster burst sources. While the 



inferred (intrinsic) isotropy of the burst radiation (section 
§3.3) allows us to at least eliminate that contribution to 
uncertainties in distance estimates, the additional system- 
atic error contributed by the observed scatter in the peak 
burst fluxes is still smaller than the usual other uncer- 
tainties due to the unknown neutron star mass and at- 
mospheric composition. Furthermore, without a detailed 
understanding of the degree of reprocessing occurring in 
the region around the neutron star, we cannot at this time 
determine the intrinsic peak luminosity of the radius ex- 
pansion bursts, from the broad distribution we have ob- 
served. 

Nevertheless, we now calculate a probable range for the 
distance to 4U 1728—34, given plausible values for the neu- 
tron star mass and atmospheric composition. We identify 
the minimum peak flux of the radius expansion bursts as 
the best estimate of the Eddington limit; this burst will 
have the smallest contribution due to reprocessed radia- 
tion, and thus will provide the best estimate of the intrin- 
sic maximum flux. Since the peak flux is typically reached 
near the end of the radius contraction, we calculate the 
gravitational redshift parameter at the neutron star radius 
-Rns = 10 km. We also reduce our observed fluxes by 20% 
to correct for the observed systematic flux offset measured 
for RXTE (see §2.2), so that the inferred Eddington flux 
is 6.2 x 10~ 8 erg cm -2 s -1 . Thus, for a 1.4(2.0) M Q neu- 
tron star with cosmic atmospheric abundance (X — 0.7), 
the distance is 4.4(4.8) kpc. For a pure He atmosphere the 
distance is 30% greater. These values are roughly consis- 
tent with previous estimates (van Paradijs 1978; Basinska 
et al. 1984; Kaminker et al. 1989) and place the source 
within 12 pc of the Galactic plane, about 4 kpc from the 
center. 
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